Dynamical Interpretation of the Nucleon-Nucleon 
Interaction and Exchange Currents in the Large Nc Limit 

D.O. Riska* 

Helsinki Institute of Physics and Department of Physics 
00014 University of Helsinki, Finland 



Abstract 

Expression of the nucleon-nucleon interaction to order 1 /Nc in terms of Fermi Invariants 
allows a dynamical interpretation of the interaction and leads to a consistent construction of 
the associated interaction currents to order I/Nq- The numerically significant components 
of 4 different modern realistic phenomenological interaction models are shown to admit very 
similar meson exchange interpretations in the large Nq limit. Moreover the ratio of the 
volume integrals of the leading, next-to-leading and next-to-next leading order terms in 
these interaction models is roughly 300:5-10:0.1, which corresponds fairly well to the ratios 
of l/N^ between the terms that would be suggested by the l/Nc expansion if Nc = 3. The 
Nc dependence of the electromagnetic and axial interaction currents that are associated with 
these interaction components is derived and compared to that of the corresponding single 
nucleon currents. 



* riska@pcu . helsinki . fi 



1 Introduction 



The most promising QC-D— based perspective on nuclear structure may very well be that pro- 
vided by the large color number {Nc) limit, which allows a systematic series expansion in 1/Nc 
for hadronic observables, where the first few terms appear to capture the key phenomenological 
features of the structure of the baryons The leading terms of the 1/Nc series expansion of 
the components of the nucleon-nucleon have been shown to correspond well with the strongly 
coupled boson exchange terms in a phenomenological boson exchange model for the interaction 
0. While the large Nc limit itself cannot determine the radial behavior of the potential compo- 
nents, this result nevertheless provides a systematic, if approximate connection between nuclear 
phenomenology and QCD. 

Here it will be shown that in fact all the most commonly employed modern realistic phe- 
nomenological interaction models a) allow an interpretation in terms of phenomenological boson 
exchange and b) that the numerically significant (or well determined) components of these inter- 
actions correspond very well to the leading terms in the I/Nq expansion of the nucleon-nucleon 
interaction components. The next-to-leading order terms are are shown to be numerically very 
small in comparison to the terms of leading order. The ratio of the volume integrals of the lead- 
ing, next-to-leading and next-to-next leading order terms in these interaction models is roughly 
300:5-10:0.1, which corresponds fairly well to the ratios of 1/iV^ between the terms that would be 
suggested by the l/Nc expansion if Nc = 3. 

Finally the corresponding 1/Nc scaling factors of the associated two-nucleon exchange (or 
interaction) currents are derived. The axial exchange charge operator that is associated with long 
range pion exchange is shown to scale as Nq, as the single nucleon axial charge, whereas those 
that are associated with the Fermi invariant components of the interaction scale as 1/Nc or higher 
powers in 1/Nc- The electromagnetic exchange currents scale with two additional powers of 1/Nc 
in comparison with the single nucleon current. The axial exchange current scales as 1/N^, and is 
thus smaller by 1/N^ than the axial current of the nucleon. These scaling factors correspond well 
with established nuclear phenomenology. 

The key to the dynamical interpretation of a phenomenological interaction model is to express 
it in terms of the 5 Fermi invariants S,V,T, A, P. These invariants form a unique set, that is 
independent of particle momenta, and allow an interpretation of the interaction in terms of linear 
combinations of scalar, vector, axial vector and pseudoscalar exchange mechanisms, without the 
need for explicit specification of the radial behavior of the potential functions. Moreover all of 
these invariants define unique and consistent interaction current operators for which 1/Nc 

expansions follow from those of the corresponding potentials, as will be shown below. 

The utility of the large Nc limit of QCD in nuclear phenomenology has emerged only gradually, 
the first indication being the a priori surprising phenomenological quality of the description of 
selected observables of nucleons, hyperons and nuclei in terms of Skyrme's topological soliton 
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model and its generalizations. These models represent one realization of the large A''^; limit, in 
which the baryons appear as topologically stable solutions of nonlinear chiral meson field theories 
0. In particular it was found that Skyrme's product ansatz for the baryon number 2 system led 
to a phenomenologically satisfactory description of the longest range components of the isospin 
dependent part of the nucleon-nucleon interaction 0, ^. This result is in fact very natural as these 
components of the interaction are of leading order in the 1/Nc expansion for which the Skyrme 
model should give generically adequate results. 

The modern phenomenological nucleon-nucleon interaction models that are considered here are 
the yi8 H, the CD -Bonn |10|, the Nijmegen{93) ^ and the Paris potentials. All of these 



contain a long range pion exchange tail, and have mostly or purely phenomenological short and 
intermediate range terms. Their components in the S, V, T, A, P representation are nevertheless 
remarkably similar in strength, and if parametrized in terms of single meson exchange give rise to 
"effective" meson-nucleon coupling strengths, which also are very similar. 

In the large Nc limit the transformation between the phenomenological potential components 
given in the spin representation and the representation in terms of Fermi invariants reveals a new 
result in that leading order parts of the isospin dependent scalar and vector potentials have to be 
equal in magnitude with opposite sign. If interpreted in terms of meson exchange this suggests 
that Oq and p meson exchange have equal strength in the large Nc limit. 

In section 2 the expressions for the Fermi invariant decomposition for the large Nc limit of 
the nucleon-nucleon interaction are derived. In section 3 the corresponding components for the 
4 interaction models considered are calculated. All of these are then shown to admit a meson 
exchange interaction with very similar coupling strengths. Their well determined components are 
then shown to be consistent with the interaction form suggested by the leading terms in the 1/Nc 
expansion. In sections 4 and 5 the Nc scaling factors of the corresponding electromagnetic and 
axial interaction current operators are derived. A concluding discussion completes the paper. 

2 The Nucleon-Nucleon Interaction in terms of Fermi In- 
variants in the large Nc Limit 

2.1 Phenomenological interactions 

The nucleon-nucleon interaction is commonly expressed in terms of the following set of Galilean 
invariant spin- and isospin operators: 

VNN = j:[vl + vTr'-A^„ (1) 
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where the coefficients are scalar functions and the spin operators f2j are defined as 

Qc = 1, ^Ls = L ■ S, Qt = Si2, 

nss = B^-a\ nLS2 = \{S^-L,a^-L}+ . (2) 

Because the httle group of transformations, which leave the square of the 4-velocity of the two- 
particle system invariant is the same for Galilean and Poincare transformations, the interaction 
(0) is invariant under Poincare transformations. 

The power of the leading term in the 1/Nc expansion for these potential components vf^ and 
therefore their order in 1/Nc, has been shown to be the following 0: 

Wc^c, %f^T, Vss^ss ~ 0{Nc), 

Vc^C, vf^LS, V^^T, V^S^SS, Vls2^LS2 ~ 0{l/Nc), 

vts2^LS2 ~ 0{l/N'c)- (3) 

A simple derivation of these scaling relations is given below. The numerically significant, or 
"large", potential components should therefore, on the basis of their order in 1/Nc, be the isospin 
independent central, and isospin dependent spin-spin and tensor components of the potential. This 
corresponds completely to established nuclear structure phenomenology, as the first of these terms 
corresponds to the short range repulsive and intermediate range attractive force components, and 
the latter two correspond to components of the long range pion exchange interaction. These are 
also the force components that are readily derived in the Skyrme model J^. 

While the 5 linearly independent potential invariants Clj are sufficient to describe the two- 
nucleon amplitude for nucleons on their mass shell, they are not dynamically transparent. For a 
dynamical interpretation, which allows a connection to field theory, reexpression of the interaction 
in terms of the 5 Fermi invariants S, V, T, A, P is required. These are defined as ||T^ 



^ = 1, V = -fl^l, T = ^crl^crl^, 



A = t'yhl ^ihl P = ihl (4) 

The nucleon-nucleon interaction, when expressed in terms of the Fermi invariants, takes the form 

Vnn = u{Pi)uip'2) K + ^7 ■ "(Pi)"(P2), (5) 

where the Fj, j = 1...5 represent the Fermi invariants in the order S, V, T, A, P. 

To find the relation between the potential functions Vj and the nonrelativistic potential com- 
ponents Vj in (|l]) the spin representation the Fermi invariants is required: 
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^LS2 = ■ na"^ ■ n. (6) 

Here the momenta fc, P and n are defined as k = p' — p, P = {l/2){p' + p) and n = k x P, 
where p' and p are the final and initial relative momenta. The linear relation between the spin 
operators Qj and the Fermi invariants Fj is given (to order 1/m^) in ref. |Q. Since the nucleon 
mass rriN ^ C>{Nc), only the leading terms in l/mj^ are needed for the present large Nc limit 
considerations. 

The momentum space representation of the interaction (|I]) is 

5 

Vnn = + ^7^' ■ ^i- (7) 

i 

Note that wls2^ls2 only represents part of the complete Fourier transform of vls2^ls2, which 
also contains a linear combination of a term with the operators ffi ■ P a2 ■ P — P ^c?i ■ cx2 , ^3 and 

inil . These additional terms are of the same order in 1/Nc as the term wf^g2^LS2- 

For a local potential one has explicitly 0] 

wf;{k) = An drr^jo{kr)v^{r), 
Jo 

^Lsi^) = — r / drr^ji{kr)v^s{r), 



k Jo 

^r(^) = IT drr^j2ikr)v^{r), 
k'^ Jo 

Wssi^) = 47r / drr^jo{kr)v$s{r), 
Jo 

471" /"oo 

«^L52(^) = ~J^ drr^j2{kr)vts2{r)- (8) 
2.2 Order in Nc of the potential components 

The Nc dependence of the components of the potential may be inferred directly by quark 
model considerations. Consider the single quark operators 1,(7^, r^., ajTk- Matrix elements of the 
sum over Nq such quark operators in nucleon states then depend on Nq as follows |T5[: 

Nc Nc 

{N\Y.l'\N)-Nc^ {N\Y.cr',\N) N'c 

q=l q=l 

Nc Nc 

{N\Y.rl\N) (iV|E^KI^)~^c. (9) 

q=l q=l 
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In the large Nc limit the baryon-baryon interaction may be interpreted as meson exchange, since 
the gluon lines may be replaced by qq lines in all surviving (planar gluon) diagrams. The meson- 
baryon couplings are proportional to the quark operator matrix elements above and inversely 
proportional to the meson decay constants /m, which scale like \fNc Application of the 
scaling relations @) and multiplication by to the two nucleon system directly implies that 

w^Vtc, w^VLt, Wgs^ss ~ 0{Nc), 

wc^c, w+^It, w^s^ss ~ 0(1/ Nc). (10) 

The order in Nc of the spin-orbit interaction may be found as follows. The non-local momentum 
operator P in the spin-orbit interaction operator Qls always appears in the combination P/rriM, 
where rriN ~ Nc- The isospin independent term w^^flis contains the spin-operator of one nucleon 
in combination with P/mi^, and therefore the vertex scales as 1/Nc- If the quark coupling at 
the other nucleon line is ~ 1'', it follows that wj^gQis ~ 0{1/Nc) once the overall factor 1/ 
is taken into account. For the isospin dependent spin-orbit interaction the scaling factor or the 
vertex with a spin-operator multiplying P/rriN and an isospin operator is N^. In this case the 
vertex at the other nucleon line also contains an isospin factor, so it also scales as N^. Thus, it 
follows after multiplication with that wJ^gflLS ~ 0{1/Nc)- 

To find the Nc scaling of the isospin independent quadratic spin-orbit interaction one notes that 
there is a combination of a spin operator with P/m^v at both nucleon lines, and thus the scaling 
factor is l/N^,, which after multiplication with 1/ flf gives rise to the scaling w~^g2^LS2 ~ ^i^/^c)- 
The corresponding scaling factor for the isospin independent quadratic spin-orbit interaction is 
wlgQLS2 ~ 0{^/Nc), because the isospin operators at both nucleon lines bring an additional 
factor Nc each. This completes the derivation of the Manohar-Kaplan scaling relations (|^) 0. 

The linear relation between the Fermi invariant potential components and the components 
wf is given in ref. 0]. In the present application, only the leading terms in (i.e. in 1/N^) 

have to be retained in the transformation matrix. 

To leading order in 1/Nc the Fermi invariant potential components vf are then given as: 

+ 3 , 

VZ = -Wn W 
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C 'TT'^LSi 



% = ^Wls + —Wt + -wss + -^Wls2, 

m^^ k"^ 1 m 



N - _ -L _ iibj^rv 
Wro W^r WoQ Wi 



1,2 1.2 1.2 1 2 1.2 

^ + ^ + +1+1 '"'N'^ + 

32m% ^ 16 2 ^ 2 4 ^^2' 
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fc^ _ 1 _ m%k'^ 
% = y^T + 2^ss + -^—^LS2, 

1.2 1.2 1.2 1 2 1.2 

f I = ^r'W^r H ■w^f'^ H 'i^T H — w^Ts : — ""^7^, 

^ 32m% ^ 16 2 ^ 2 4 

A;^ _ 1 _ ■m%k'^ 

= y^T + 2^55 ^^L52, 

-, 9 

m 



fp = VZmljW^ . (11) 

From these expressions it follows that both the isospin independent and the isospin dependent 
scalar and vector potentials vf S, V are ~ 0{Nc). The isospin dependent tensor and axial vec- 
tor potentials T, A are ~ 0{Nc), whereas the corresponding isospin independent interaction 
components T, A are~ 0{1/Nc)- The orders of the pseudoscalar potentials are VpF^^ Nq 
and Vp P ^ N^. All the coefficients in the expression of the pseudoscalar operator P in terms of 
the spin operators Qj are ~ 1/'"^^? @]- The corresponding potentials in the spin representation are 
therefore smaller by l/N^; and as a consequence the highest order in Nq of the potentials compo- 
nents in the spin representation is 0{Nc)- This becomes explicit if the pseudoscalar invariant P 
is replaced by the on-shell equivalent pseudovector invariant P': 

P' = ^llKll iX^l (12) 

which has the form required by spontaneous breaking of the approximate chiral symmetry of 
QCD. Because the pseudovector invariant builds in the required pair suppression in the pion- 
nucleon amplitude, and vanishes with the 4-momentum of the exchanged system, it is simpler 
to apply in relativistic approaches to nuclear reaction phenomenology |1^. These results are 
summarized in Table 1. 



2.3 Dynamical interpretation of the potential components 

The components of highest order in 1/Nc of the Fermi invariant potential components in Table 1 
admit a direct dynamical interpretation in terms of meson exchange between the nucleons. 

The only isospin independent potential components that are of order Nc (i.e. "large") are the 
scalar and vector components and Vy as well as the pseudoscalar component Vp. The first one of 
these corresponds to the largest component in the two-pion exchange interaction between nucleons 



17, 0, which commonly is modeled in terms of a strong scalar ("a") meson exchange mechanism 



19, 20|. This potential component gives rise to an attractive interaction at intermediate range 



and the main agent for nuclear binding. The second corresponds to the short range repulsion 
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Isospin 


S 


V 


T 


A 


P 


1 


Nc 


Nc 


l/Nc 


l/Nc 


Nc 




Nc 


Nc 


Nc 


Nc 


Nh 



Table 1: Order of the leading term in the l/Nc expansion of the Fermi invariant potential com- 
ponents. 



between nucleons, and is commonly modeled in terms of an ui meson exchange interaction, with 
an overstrength effective u nucleon coupling constant ||20[. The last (pseudoscalar) term, which 



may be interpreted as r] meson exchange is large in the earlier potential models, but very small in 
the most recently developed models. 

Among the isospin dependent Fermi invariant potential components in Table 1 the pseudoscalar 
component is of order N^ and are thus the largest of all terms in the l/Nc counting scheme. This 
strong pseudoscalar exchange component Vp is immediately interpretable in terms of the long 
range pion exchange interaction between nucleons, which is expectedly "strong" , more because of 
its long range than because of the strength of the pion-nucleon coupling. 

That the other isospin dependent Fermi invariant potential components are of order Nc also 
corresponds to established nucleon-nucleon phenomenology. To see this, it is worth noting that 
an isospin 1 vector meson (p) exchange interaction may be expressed in terms of Fermi invariants 



as follows 



k{1 + k)P ^ ^ k{1 + k) plp2p 

Here to leading order in l/m% (or l/N"^): 



(13) 



PlPl = -{ml + ^ + 2P^). (14) 
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Phenomenological boson exchange interaction models typically contain a p meson exchange in- 



teraction, where the "effective" tensor coupling of the p meson to nuclei is large k ~ 7 even 
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though the vector couphng is smalL Because of the large tensor coupling such an isospin dependent 
vector meson exchange interaction contributes strongly to all the Fermi invariant potential com- 
ponents, except the axial vector invariant. The order Nc of the potential components w^, fy, 
thus may be interpreted in terms of a strong p meson exchange (or p exchange-like) interaction. 
The isospin independent vector meson {uj) exchange in contrast only contributes strongly to the 
vector V invariant potential, because of the small u;— nucleon tensor coupling. 

While the isospin dependent pseudoscalar potential f p P is of order the corresponding 
isospin independent pseudoscalar exchange potential fp P is only of order Nc- This corresponds 
to the phenomenological finding in boson exchange models for the nucleon-nucleon interaction 
that the rj meson exchange term is mostly much weaker than the tt meson exchange interaction 
component. 



3 Large Nc components of phenomenological interaction 
models 

3.1 The scalar potential 

For the visualization of the Fermi invariant potential components it is useful to plot the poten- 
tial components in configuration space. The expression for the isospin independent and isospin 
dependent scalar potentials in configuration space are obtained from (0) as 

4(0 = 4^^(0 + ^/ dr'r'vtsir'l 

vsir) = dr'r'vZsir') - (r) - 3 ^ dr' ^} + - 5^,(r) 



rrir 



oo 



-^{r'vls2ir) - 3 dr'r' vls2ir')}, (15) 

Here, and below, only terms of leading order in 1/Nc have been retained, so that all terms 
on the r.h.s. are of the same order in 1/Nc- These potential components as obtained from 
the phenomenological ^18 |p, the CD — Bonn the Nijmegen{93) [|ri[] and the Paris |jl2| 
potentials are shown in Figs. 1 and 2. 

In the case of the V^18 potential, the operator form of which contains a quadratic spin-orbit 
interaction of the form (L ■ S")^, we employ the relation 

(L ■ Sf = 2P + 2fii52 - 2hLs (16) 

to reduce the interaction to standard form (|1]) in combination with terms of the form L^. As the 
order of the remaining terms L^, (r^ ■ t^)L'^, {ff^ ■ (?^)L^ and (r^ ■ t^){&^ ■ (t'^)L'^ in the V18 potential 



9 



may, by the argument in section (2.2), be shown to be at most of order and 
respectively (once multiphed by we shall not have to retain them in the leading 

order analysis here. The spin-orbit {^Ils) and quadratic spin-orbit {flLS2) interaction terms in the 
1^18 interaction shall however be included with the corresponding terms in the expressions (pIS]). 

The CD — Bonn interaction model is a single boson exchange model that is expressed in terms 
of Fermi invariants and in the case of vector meson exchange in the form, which is reducible to 
Fermi invariants by eqn. (p^Sf ). Both the Nijmegen{93) and the (parametrized) Paris interaction 
models are expressed in the standard form (|l]). The central components of the Paris potential 
do in addition contain velocity dependent terms with the operator form P^. As these, by the 
argument in section (2.2) are down by two orders in 1/Nc in comparison to the local terms, they 
are not considered here. 

In Table 2 are listed the volume integrals of the the leading terms in the 1/Nc expansion of 
the Fermi invariant potential components of the phenomenological interaction models. In the case 
of the isospin independent scalar interaction these volume integrals are all of the order 10 fm^, 
ranging from -6.3 fm^ for the Paris potential to -13.5 for the CD — Bonn potential. 

If the isospin independent scalar potential component is interpreted as due to a single scalar 
meson exchange interaction, the volume integral would equal —gg/m'g, where gs is the scalar 
meson - nucleon coupling constant, and ms the scalar meson mass. The mass of the lowest scalar 
meson (a or /o(400 — 1200)) is not well established [0, but it is typically taken to be of the order 
550 - 750 MeV in phenomenological nucleon- nucleon interaction models . If is taken to 
be 600 MeV, within this range, the value of the effective scalar meson coupling constant for these 
interactions range between 7.6 {Paris) and 11.2 {CD — Bonn) (Table 3). These "effective" scalar 
meson coupling constant values should be viewed as lower limits, as any short range form factors 
multiplying the scalar meson exchange interaction would increase these values. It is striking how 
closely similar values obtain for the "effective" scalar meson coupling constants with the 4 different 
interaction models, two of which do not contain any scalar meson like term in their parametrized 
forms. 

The isospin independent scalar interaction components of the 4 interaction models are rather 
similar for nucleon separations beyond 0.5 fm. The radial shapes differ considerably at short dis- 
tances, however, ranging from attractive to repulsive. The two boson exchange models {Nijmegen{93) 
and CD — Bonn are repulsive by construction. The more phenomenological V^18 and Paris in- 
teraction models are repulsive at short distance, a feature which cannot be interpreted in terms 
of meson exchange, and which otherwise only is known from the Skyrme model p. 

The isospin dependent scalar interaction components of the 4 phenomenological interaction 
models are shown in Fig. 2. These interaction components are of order 1/Nc and should therefore 
be weaker by and order of magnitude than the isospin independent scalar interaction components. 
This is in fact also revealed by comparison of Figs. 1 and 2. The volume integrals of these 
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interaction components range between -3.2 and -4.4 fm^, as shown in Table 2, and are thus about 
a factor 3-4 smaller than those of the corresponding isospin independent interactions. 

In a meson exchange interpretation the isospin dependent scalar interaction arises from ex- 
change of scalar mesons with isospin 1. The lightest of these is the ao(980) In a simple meson 
exchange model the volume integrals of these interactions should equal —gi^/m^ao^ where Qa^ is 
the ao(980)— nucleon coupling constant and is the mass of the ao(980). If the value of the 
"effective" ao(980) nucleon coupling constant is extracted from these volume integrals, the values 
range from Qa^^ = 9.0 to Qa^ = 10.4 (Table 3). In this case all the interaction models give very 
similar values for this coupling constant. 



3.2 The vector potential 

The expression for the isospin independent and isospin dependent vector potentials in configuration 
space are obtained from (|TT]) as 

^y(^) = ^^cW - dr'r'vtsir'), 



^y(^) = — — 1^ dr'r'vLsir') + -{vTir) -3 j dr' —^} - -Vss{r) 



m 



+ ^{r'vls,{r) - 3 dr'r' vls,{r')}, (17) 
These potential components as obtained from the phenomenological V^18 0, the CD — Bonn [jlO 



the Nijmegen{93) [|n[] and the Paris potentials are shown in Figs. 3 and 4. Note that in 
the large Nc limit Vy{r) = —Vg{r). In a boson exchange interpretation this implies equality in 
strength between the ao(980) and p(770) exchange interactions. 

The isospin independent vector interaction admits an obvious interpretation in terms of a 
strongly repulsive u meson exchange interaction. The shape of this interaction component is very 
similar for all the interactions considered, and the volume integrals are also remarkably similar, 
ranging from 8.7 {Nijmegen{93)) to 11.6 fm^ [CD — Bonn). If the "effective" tu— nucleon coupling 
constant is extracted from the volume integrals, by setting them equal to g'^/m'^ the values range 
between 11.7 and 13.5. The similarity between these coupling constant values is particularly 
notable in the case of the ^18 and the parametrized Paris interaction models, as neither one of 
these contains any explicit vector meson exchange term. 

The isospin dependent vector interaction, by eqn. ([T3|) , admits an interpretation in terms of 
p— meson exchange. This interaction component is roughly 3 times weaker than the corresponding 
isospin independent interaction. The volume integrals of the isospin dependent vector interaction 
range from 2.9 to 4.4 fm^ for the 4 phenomenological interactions considered here. In a simple 
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Component 


V18 


CD - Bonn 


Nijmegen{93) 


Paris 


0{Nc) 


vt 


-8.7 


-13.5 


-10.3 


-6.3 


Nc 




-3.2 


-3.3 


-3.3 


-4.4 


Nc 




9.4 


11.6 


8.7 


10.2 


Nc 




3.2 


3.2 


2.9 


4.4 


Nc 




-0.1 


0.0 


0.0001 


0.03 


l/Nc 




0.6 


0.1 


0.001 


0.46 


Nc 




-0.1 








0.03 


l/Nc 




0.6 








0.46 


Nc 


vt 


9.8 


0.0 


0.35 


18.0 


Nc 


Vp 


360 


338 


323 


352 


Nh 



Table 2: Volume integrals (in fm^) of the leading terms in the l/Nc expansion of the Fermi 
invariant potential components phenomenological interaction models. 
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Component 


V18 


CD - Bonn 


Nijmegen{93) 


Paris 


9a 


9.0 


11.2 


9.8 


7.6 


9ao 


9.0 


9.0 


9.0 


10.4 


9u> 


12.2 


13.5 


11.7 


12.7 




7.0 


7.0 


6.3 


10.1 


9n 


8.7 


0.0 


1.8 


11.7 


9n 


13.4 


13.0 


12.7 


13.2 



Table 3: Effective meson-nucleon coupling values tfiat correspond to pfienomenological interaction 
models. 
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p— meson exchange model these volume integrals have the form (1 + K)(yf^/m^, where Qp and k are 
the vector and tensor p— meson-nucleon coupling constants. The canonical value for the p— nucleon 
vector coupling constant is gp/4:7r ~ 0.5. If this value for Qp is used in the expression for the volume 
integral, an "effective" value for Kp can be determined. The volume integrals in Table 2 for the 4 
interaction models considered then give values for Kp in the range 6.3 - 10.1. These values support 
the early finding that Hp ~ 6.6 ||2^. Both the Nijmegen93 and the CD — Bonn interaction use 
the somewhat larger value gp/4:TT = 0.84 for the p-nucleon coupling constant. This larger value 
would reduce the values for Kp in Table 3 by about 40 %. 

From eqns. (|15D and (0) it follows that in the limit of large Nc the isospin dependent scalar 
and vector exchange potentials should be equal in strength and of opposite sign. If imposed on a 
boson exchange model this condition implies that 

gl = (1 + Kp)gl. (18) 

Comparison of the form and strength of these potential components for the CD — Bonn and 
Nijmegen{93) potential models in Figs. 2 and 4 shows that these potential models satisfy this 
constraint very well. 



3.3 The tensor potential 

The expression for the isospin independent and isospin dependent tensor potentials in configuration 
space are obtained from (p!T| ) as 



1 ;/ 2 / 1 / 

SZmif r 16 

\ roo V~(t'^ 1 

v¥{r) = --{d^ir)-3j^ dr' + -dssir) 



m 



oo 

{r' vls2{r) - 3 dr' r' vls^^r')} , (19) 



These potential components as obtained from the phenomenological V^18 the CD — Bonn |10 
the Nijmegen{93) |TI[] and the Paris [|I^] potentials are shown in Figs. 5 and 6. 

In a meson exchange model the tensor potential arises from vector meson exchange when 
the vector mesons couple to the nucleon with a Pauli (tensor) coupling. As in the CD — Bonn 
interaction model the isospin independent u meson exchange interaction has no Pauli coupling 
term Vj^ vanishes in this interaction model. 
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By the Nc counting rules in Table 1 the isospin independent tensor interaction should be 
smaller by l/N"^ than the isospin dependent tensor interaction. Comparison of the corresponding 
set of volume integrals in Table 2 shows that the phenomenological interaction models considered 
here satisfy this rule well. 

The order Nc of the isospin dependent tensor interaction alone does not explain why this 
interaction is one order of magnitude smaller than the corresponding vector interactions for all 
the phenomenological potential models. For boson exchange models the reason is readily seen in 
eqn. (|13]), which shows that the tensor coupled vector meson exchange interaction is suppressed 
by an overall factor which is only partially counteracted by the large tensor coupling Kp. 

For the phenomenological interaction models the reason for the weakness of the isospin dependent 
tensor interaction is to be found in the fact that pion exchange gives rise to the bulk of the isospin 
dependent spin-spin and vector interactions, and this pion (or pseudoscalar) exchange contribution 
is exactly cancelled in the combination of v^g and in eqn. ([T9|) . This latter argument may also 
be extended to the isospin independent tensor interaction, where the pseudoscalar - in this case 
the f] meson - exchange contribution cancels in the combination of v~gg and in ([1^) . 



3.4 The axial vector potential 

The expression for the isospin independent and isospin dependent axial vector potentials in con- 
figuration space are obtained from ([Tl|) as 



1 // 2 / 1 / 



,2 



1 /"OO 7? I I 1 Tn 

W = -^i^rir) - dr' + -Vssir) 



2 

mi 



+ ^{r'vls2ir) - 3 / drV vls2ir')}, (20) 



These potential components as obtained from the phenomenological V18 and the Paris |T2 
potentials are shown in Figs. 7 and 8. In a single meson exchange model only axial vector meson 
exchange gives rise to an axial vector interaction. Since the CD — Bonn and the Nijmegen93 
boson exchange interaction models do not contain any oi meson exchange terms these interaction 
models do not have any axial vector exchange components. The phenomenological V"18 and Paris 
potential models give rise to small A interaction components. These satisfy the Nc counting rules 
in Table 1, by which the isospin dependent axial vector interaction should be larger by Nq than 
the isospin independent one as may be seen from the volume integrals in Table 2. The substantial 
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strength of the axial vector exchange interaction components of the Paris potential model are 
likely to be artifacts of the parametrization rather than genuine consequences of the NN scattering 
data. 



3.5 The pseudoscalar potential 

The expression for the isospin independent and isospin dependent pseudoscalar potentials in con- 
figuration space are obtained from ([TT|) as 



v+{r) = -v+{r) - ^ / dr'r'vtsir') + 12m^ / dr'r' / dr" -^^-^ 



4 . 2 



r Jr 



, oo poo 7) ( 



Vp{r) = Vlm% / dr'r' / dr" ' \^ \ (21) 

Jr Jr' r 

These potential components as obtained from the phenomenological 1^18 0, the CD — Bonn and 
the Nijmegen{93) |ll|] and the Paris [|1^] potentials are shown in Figs. 9 and 10. 



In a single meson exchange interpretation the longest range mechanism, that contributes to the 
isospin independent pseudoscalar interaction, is rj meson exchange. This interaction component is 
not well constrained by nucleon-nucleon scattering data. As a consequence the phenomenological 
interaction models considered here give widely different results for this interaction component, as 
is evident in Fig. 9. 

The 77— nucleon coupling constant is also not well known. An estimate for this coupling constant 
may be obtained from the volume integrals of the phenomenological interaction models if these 
are set to equal g'^/rn^. From the volume integral values listed in Table 2 one then obtains values 
for 51^, which range from 1.8 {Nijmegen{93)) to 11.7 (Paris) (Table 3). Analyses of observables, 
other than nucleon-nucleon scattering, as e.g. r^— meson photoproduction suggest that the coupling 
constant value should not exceed 2.2 p3| . 

The isospin dependent pseudoscalar exchange interaction is strong and has long range. Its main 
component is the long range pion exchange interaction, which is built into all the phenomenological 
interaction models considered here. Because of this all the interaction models converge for nucleon 
separations larger than 1 fm. This interaction component is also the strongest component by Nq 
counting, as it scales as A*"^ (Table 1). The volume integrals of this interaction component are 
listed in Table 2 for the 4 interaction models considered, and range from 323 fm^ to 360 fm^. The 
corresponding values for the pseudoscalar pion-nucleon coupling constant g.^^ range from 12.7 to 
13.4, when extracted by equating the numerically determined volume integrals with gl/rn^ (Table 
3). 
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4 Electromagnetic exchange current in the large Nc limit 



The isoscalar and isovector components of the single nucleon charge operator are of order Nq 
and A^^ respectively, as may be inferred from the scahng relations (P). Because of their inverse 
dependence on the nucleon mass the isoscalar and isovector convection current components of the 
single nucleon current scale as A^^ and A'",^^ in comparison. In the case of the spin component 
of the single nucleon current the isoscalar component scales in contrast as A*",^^, whereas the Nc 
dependence of the isovector part is N^. 

Associated with the isospin dependent interactions expressed in terms of Fermi invariants 
are corresponding two-nucleon electromagnetic interaction currents with the isospin dependence 
(r^ X r^)3 1^. In a meson exchange interpretation, these interaction currents correspond to the 
currents carried by exchanged charged mesons (Fig. lib). These current operators are listed in 
Table 4, along with the corresponding Nc scaling power. The Nc scaling factors of these current 
operators may be derived by reduction to the spin representation, and application of the relations 
In the expressions in Table 4 the potential functions vf{k) represent the Fermi invariant 
potential components in momentum space. The Nc dependence of these interaction currents 
should be compared to the Nc scaling of the corresponding components of the single nucleon 
charge and current operators. 

The charge components of the interaction currents in Table 4 are suppressed by their de- 
pendence on the energy exchanged between the nucleons. This is proportional to the nucleon 
momenta and inversely proportional to the nucleon mass (~ 1/Nc)- In a frame with no energy 
transfer they vanish. In the reduction to the spin representation exchange charge operators also 
arise in the elimination of negative energy intermediate states ("pair terms"). These are listed in 
Table 5. These operators have an overall factor l/m%, which arises from the small components of 
the Dirac spinors and the propagator of the intermediate negative energy propagator (Fig. 11a). 
The order in 1/A^c of these operators follow from the general scaling rules (P) and the overall 
factor 1/m^ ~ V^c- Table 5 and in subsequent tables, the momenta of the two nucleons are 
denoted by Pj = l/2{pj' + pj) with j = 1,2 respectively. 

The largest Nc scaling factors of the exchange charge operators that are associated with the 
5 Fermi invariant components of the interaction are also listed in Table 5, separately for the 
isospin independent and isospin dependent interaction components. The order of the exchange 
charge operators are smaller by I/Nq in comparison with the corresponding single nucleon charge 
operators. The terms with the highest Nc dependence are those, which are independent of spin and 
isospin and those, which contain the spin-isospin bilinears aiTj of both nucleons. The remaining 
terms are smaller by at least one power of 1/Nc- 

The interaction currents that are listed in Table 4 arise from the coupling of the electromagnetic 
field to the charge exchanged between the nucleons (Fig. lib). Their explicit expressions in the 
spin representation have been given in ref. 0]. The components of these currents, which are of 
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Fermi invariant 






a 
J 
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IH^ X ^ J3 kl-kl 

+ 7^(7^ -A^i)- 7^(7^ -M]} 


Nc' 


T 


[{k2^, - kif,)T + k2u(yl^al^ - ki^aj^^a^} 


Nh 


A 


^{r' X r%'-^^^p^{k2, k,,)A 


Nh 


P 




Nc' 



Table 4: Interaction currents and Nc scaling factors associated with the interaction components. 
The fractions of the imparted momentum q to the two nucleons are denoted ki and Note that 
in the case of the A invariant a conserved current obtains only in combination with a P invariant 
term Q]. 
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Fermi invariant 
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[vi{k2)ei + v2{k2)i2] + {I ^ 2) 
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[vUk2)ei + Vp{k2)e2] + {I ^ 2) 







Table 5: Exchange charge operators and the highest Nc scahng factors associated with the inter- 
action components. The fractions of the imparted momentum q to the two nucleons are denoted 
ki and k2 respectively. The isospin factors are defined as ei = (1 + T3 )/2 and 62 = (r^ • + r|)/2. 
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highest order in Nc are hsted in Table 6. 

The Nc dependence of these internal coupling interaction currents that are associated with 
the different Fermi invariant components of the phenomenological interaction model is smaller by 
one power of 1/Nc than that of the isovector part of the single nucleon current, with exception for 
the currents that are associated with the T and A invariants, which anomalously scale as A^"^. As 
these interaction components are very weak in all the phenomenological models considered, and 
are suppressed by l/m^ ~ l/^c- meson exchange models, it is natural to conclude that their 
Nc dependence in practice also is iVj^ rather than iV^. 

In the reduction to the spin representation the negative energy pole terms give rise to seagull 
type exchange currents, in which the electromagnetic field couples to either one of the interaction 
nucleons with a point interaction to the interaction (Fig. 11a). The exchange currents of this 
type, which have the highest Nc scaling factors (in this case 1/Nc) are listed in Table 7. The 
1/Nc dependence of these current operators are similar to the internal charge coupling currents 
in Table 6. 

5 Axial exchange current in the large Nc limit 

The axial current operator of a single nucleon, 

A± = -QA^Ti, (22) 

is of order A*"^ (cf. (^). The axial vector interaction currents that are associated with the Fermi 
invariant components of the nucleon-nucleon interaction have been listed in ref. [^. These 
currents only arise in the reduction to the spin representation, which brings in an overall factor 
1/M^. The 1/Nc scaling factors of these axial exchange current operators may be determined 
from their operator structure by means of the relations (^. It follows from their vector nature 
that they are 1/Nc s^ialler than those of the corresponding interaction components. 

The subset of axial exchange currents that are associated with the Fermi invariant interaction 
components, and which are of order 1/Nc listed in Table 8. In the case of the pseudoscalar in- 
variant the reduction to the spin representation brings an overall factor 1/Mf^. The axial exchange 
current that is associated with the pseudoscalar invariant is therefore of the order N^^. 

The axial exchange charge operators that are associated with Fermi invariant decomposition 
of the nucleon-nucleon interaction have been derived in ref. and are listed in Table 9. For 
these operators the reduction to the spin representation only brings in an overall factor of l/m%. 
The 1/A"c scaling factors of these operators are accordingly smaller by a factor 1/A^^ than those 
of the corresponding interaction components. 



20 



Fermi invariant 
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Table 6: Interaction currents that are associated with the Fermi invariant interaction components 
that arise from e.m. couphng to the exchanged charge, which have the highest Nc scahng factors 
associated with the interaction components. The fractions of the imparted momentum q to the 
two nucleons are denoted ki and /c2 respectively. 
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Fermi invariant 
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Table 7: Interaction currents that arise from coupling to intermediate negative energy pole terms 
that have the largest Nc scaling factors. The fractions of the imparted momentum q to the two 
nucleons are denoted ki and k2 respectively. 
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Fermi invariant 
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Table 8: Axial exchange current operators with the highest Nc scaling factors associated with the 
interaction components. 



23 



The axial charge operator of a single nucleon 



Al = -g^^T^^ (23) 
rriN 

is readily seen to be of order A*"^ by means of the relations (P). The main axial exchange charge 
operators listed in Table 6 are in comparison of order N^^. This may be seen by referring to the 
operator scaling relations @. 

There is no axial exchange charge operator that would be directly proportional to the pseu- 
doscalar invariant potential component. There is however an axial exchange charge operator, 
which is of order N^, and which arises from the long range pion exchange interaction, which gives 
rise to the bulk of the pseudoscalar interaction component. This axial charge operator was derived 



in ref. lit], and has the form 



n A m 1 p-"i,r'"i2 

Alin) = '^i^na^ + a') ■ n,(r- x r%il + ^f- . (24) 

4vr 2/^ m^ri2 m^ri2 

Because the pion decay constant ~ Nq in the large Nq limit it is readily seen by the relations 
(y) that this operator is of order N^. 

The contributions of the two-nucleon exchange current terms to the axial charge operator are 
conventionally expressed in the from of an effective single nucleon axial charge operator. The 
exchange current contributions then represent nuclear enhancement factors of the single nucleon 
axial charge operator (^31) . Experimentally this enhancement factor is a factor 2 as measured by 
first forbidden (3 transitions of nuclei in the lead region About one half of this enhancement 
may be attributed to the pion exchange operator (0), which is natural as its order in the 1/Nc 
expansion is the same as that of the single nucleon operator ( pSj) . The bulk of the remainder arises 
from the isospin independent scalar contribution in Table 9 [Q]. While the order of this term is 
1/Nc, it gives rise to a direct term matrix element, which involves an additional factor of Nq in 
the nucleon density, and therefore its contribution to the "effective" single nucleon matrix element 
is also of order N^. 



6 Discussion 

This analysis of 4 commonly employed phenomenological nucleon-nucleon interaction models re- 
veals that the structure of their components is completely consistent with their corresponding 
dependence on Nc (or 1/Nc). For all the interaction models the isospin dependent pseudoscalar 
components, which contain the long range pion exchange interaction, have volume integrals that 
are larger than those of any other interaction component by more than an order of magnitude. 
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Fermi invariant 
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Table 9: Axial exchange charge operators and highest Nc scaling factors associated with the 
interaction components. The fractions of the imparted momentum q to the two nucleons are 
denoted ki and k2 respectively. 
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This interaction component also scales with the largest power of A^^^; (A*"^) (Table 2). As this in- 
teraction component shows very little variation between 3 of the 4 interaction models considered, 
it may be viewed as well determined, if not completely settled in form (Fig. 10). 

The interaction components that follow the isospin dependent pseudoscalar interaction in 
strength are the isospin independent scalar and vector interactions, which scale as Nc- These 
interaction components are those responsible for nuclear binding (f J) and short range repulsion 
between nucleons {vy). For distances shorter than 0.4 fm the variation in form of these interaction 
components between the 4 considered phenomenological interaction models is substantial. 

Next in strength are the isospin dependent scalar and vector interactions, which also scale as 
Nc- The volume integrals of these interaction components are about half as large as those of 
the corresponding isospin independent interactions. A new finding is that these two interaction 
components have to have equal magnitude and opposite sign in the large A'c limit. 

The tensor and axial vector interaction components are very weak for all the considered phe- 
nomenological interaction models. For the isospin independent tensor and axial vector components 
this is completely consistent with their Nc dependence, which is 1/Nc- The smallness of the isospin 
dependent tensor and axial vector exchange interaction components cannot be explained by their 
Nc dependence alone, as they scale as Nc- In meson exchange models the smallness of is 
however natural, as it arises from tensor coupled vector mesons, and contains an overall factor 
l/m% (0). 

The least well understood interaction component is the isospin independent pseudoscalar inter- 
action (Fig. 9). This component, which scales as A"^, is vanishingly small in two of the considered 
phenomenological interactions and even stronger than Vy in the other two. The longest range 
part of this interaction component arises from t] meson exchange. The large variation between 
the phenomenological interaction models for this interaction component reflects the continuing 
uncertainty concerning the strength of the r^— nucleon coupling. The ^18 and Paris interaction 
models have strong "effective" 1] nucleon couplings, which correspond to the old pseudoscalar 
SU{3) coupling model by which gnNN = g-KNN/VS. 

The present results for the A'c dependence of the exchange currents follow those of the corre- 
sponding interactions. In the case of the electromagnetic interaction current the result that they 
are smaller by 1/A^^ than the corresponding interactions fits with the phenomenological finding 
that electromagnetic exchange current contributions to nuclear observables typically represent 
effects of the order 10%, the exceptions being situations where single nucleon current matrix 



elements are suppressed |24] 



The interaction current contributions to the axial charge operator are known to be significant in 
comparison with the axial charge operators of the single nucleons This feature also emerges 
from the 1/Nc expansion, as the axial exchange charge operators scale in 1/Nc as the single 
nucleon charge operator, or with one more power of 1/Nc- 

The overall conclusion is thus that the ordering of nuclear interactions and interaction current 
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operators by their dependence on I/Nq corresponds well with their significance ordering based on 
empirical evidence and phenomenlogical analysis. 



Acknowledgments 

The hospitality of the W. K. Kellogg Radiation Laboratory of the California Institute of 
Technology is gratefully acknowledged. I am indebted to C. Helminen for a careful reading of the 
manuscript. Research supported in part by the Academy of Finland by grant 54038. 

References 

[1] E. Jenkins, Ann. Rev. Nucl. Sci. 48 (1999) 41 

[2] D. B. Kaplan and A. V. Manohar, Phys. Rev. C56 (1997) 76 

[3] M. Kirchbach, D. O. Riska and K. Tsushima, Nucl. Phys. A542 (1992) 616 

[4] K. Tsushima, D. O. Riska and P. G. Blunden, Nucl. Phys. A559 (1993) 543 

[5] T. H. R. Skyrme, Proc. Roy. Soc. A 260 (1961) 127 

[6] E. M. Nyman and D. O. Riska, Reports on Progress in Physics 53 (1990) 1137 

[7] T. H. R. Skyrme, Nucl. Phys. 31 (1962) 556 

[8] D. O. Riska and E. M. Nyman, Phys. Lett. B183 (1987) 7 

[9] R. B. Wiringa, V. G. J. Stoks and R. Schiavilla, Phys. Rev. C51 (1995) 38 
[10] R. Machleidt, Phys. Rev. C63 (2001) 024001 
[11] V. G. J. Stoks et al., Phys. Rev. C49 (1994) 2950 
[12] M. Lacombe et al., Phys. Rev. C21 (1980) 861 
[13] M. L. Goldberger et al., Phys. Rev. 120 (1960) 2250 
[14] P. G. Blunden and D. O. Riska, Nucl. Phys. A536 (1992) 697 
[15] A. Wirzba, M. Kirchbach and D. O. Riska, J. Phys. G20 (1994) 1583 

27 



[16] J. A. Tjon and S. J. Wallace, Phys. Rev. C35 (1987) 280 

[17] M. H. Partovi and E. L. Lomon, Phys. Rev. Lett. 22 (1969) 438 

[18] M. Chemtob, J. W. Durso and D. O. Riska, Nucl. Phys. B38 (1972) 141 

[19] D. O. Riska and G. E. Brown, Nucl. Phys. A53 (1970) 8 

[20] R. Machleidt, K. Holinde and Ch. Elster, Phys. Rept. 149 (1987) 1 

[21] D. E. Groom et al., European Phys. J., C15 (2000) 1 

[22] G. H5hler et al, Nucl. Phys. B144 (1976) 505 

[23] L. Tiator, C. Bennhold and S. S. Kamalov, Nucl. Phys. A580 (1994) 455 

[24] D. O. Riska, Phys. Rept. 181 (1989) 208 

[25] K. Tsushima and D. O. Riska, Nucl. Phys. A549 (1992) 313 

[26] K. Kubodera, J. Delorme and M. Rho, Phys. Rev. Lett. 40 (1978) 755 

[27] E. K. Warburton, Phys. Rev. Lett. 66 (1991) 1823 



Figure captions 

Fig. 1 Isospin independent scalar potential coefficients for the V^18 the CD — Bonn |rD|, 
the Nijmegen{93) |ll[] and the Paris [|1^] interaction models. 

Fig. 2 Isospin dependent scalar potential coefficients for the V18 [Q, the CD — Bonn []T^, the 
Nijmegen{93) |jll| and the Paris [0 interaction models. 



Fig. 3 Isospin independent vector potential coefficients for the V18 [0, the CD — Bonn |TD 
the Nijmegen{93) |ll[] and the Paris [|1^] interaction models. 

Fig. 4 Isospin dependent vector potential coefficients for the V18 [^, the CD — Bonn |T0[, the 
Nijmegen{93) fill] and the Paris [0 interaction models. 



Fig. 5 Isospin independent tensor potential coefficients for the V18 [p, the Nijmegen{93) |TT 
and the Paris [O interaction models. 
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Fig. 6 Isospin dependent tensor potential coefficients for the V^18 [Q, the CD — Bonn |TD|, the 
Nijmegen{93) |jTl]] and the Paris interaction models. 



Fig. 7 Isospin independent axial vector potential coefficients for the V18 ^ and the Paris 



12 interaction models. 



Fig. 8 Isospin dependent axial vector potential coefficients for the V"18 0, the Nijmegen{93) 
and the Paris interaction models. 

Fig. 9 Isospin dependent pseudoscalar potential coefficients for the 1^18 0, the Nijmegen{93) 
TTH and the Paris interaction models. 



Fig. 10 Isospin dependent pseudoscalar potential coefficients for the V18 p|, the CD — Bonn 



TOU , the Nijmegen(93) [|Tl|] and the Paris [|T2| interaction models 



Fig. 11 (a) Point coupling of an extremal field to one nucleon and a nucleon-nucleon interaction, 
(b) Coupling to the charge exchanged by a nucleon-nucleon interaction. 
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